Frustrated Order by Disorder: the Pyrochlore Antiferromagnet with 

Bond Disorder 
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The classical Heisenberg antiferromagnet on the pyrochlore lattice is macroscopically and contin- 
uously degenerate and the system remains disordered at all temperatures, even in the presence of 
weak dilution with nonmagnetic ions. We show that, in stark contrast, weak bond disorder lifts the 
ground state degeneracy in favour of locally collinear spin configurations. We present a proof that 
for a single tetrahedron the ground state is perfectly collinear but identify two mechanisms which 
preclude the establishment of a globally collinear state; one due to frustration and the other due to 
higher-order effects. We thus obtain a rugged energy landscape, which is necessary to account for 
the glassy phenomena found in real systems such as the pyrochlore Y2M02O7 recently reported by 
Booth et al. jl[ to contain a substantial degree of bond disorder. 
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Geometrically frustrated magnets on the kagome 
or pyrochlore lattices, are ideal model systems for the 
study of general concepts pertaining to noncrystalline 
condensed phases of matter. Behaviour analogous to 
that of liquid || , glass J|,|||| , and ice ||,|7j phases have 
already been attributed to experimental pyrochlore sys- 
tems, while the related Gadolinium Garnet antiferromag- 
net offers liquid || and re-entrant solid behaviour ||. 
Equally, theoretical models provide the opportunity to 
study an array of many body phenomena: the classi- 
cal Heisenberg antiferromagnet on the pyrochlore lattice 
retains the behaviour of an idealized liquid right to zero 
temperature [fl3 , while in the equivalent kagome IO] sys- 
tem and in the XY antiferromagnet on the pyrochlore 
lattice |l(|[l2| thermal "order by disorder" effects lead to 
the apparition of liquid crystal like phases. 

Given this rich array of phenomena one would like to 
present geometrically frustrated magnetism as a field in 
which one can study charicatures of complex many body 
systems. However, despite some successes (see for ex- 
ample [pp|Jl^|) there remain a number of marked dis- 
agreements between experiment and theory which pre- 
vent such a desirable coherent picture. For example, 
experiments on pure Y2M02O7, which is believed to be 
well represented by a Heisenberg antiferromagnet on the 
pyrochlore lattice, show an extremely well defined spin- 
glass transition p| . Similar spin glass freezing is found in 
Tb 2 Mo 2 7 §. This is in sharp contrast to the behaviour 
of the classical Heisenberg antiferromagnet on the py- 
rochlore lattice, which shows no sign of glassy behaviour 
at any temperature |]Io|| . In fact Y2M02O7 behaves much 
more like a conventional randomly disordered spin glass 
or even the disordered compound CsNiCrFg - which has 
two kinds of magnetic ion (Cr 3+ and Ni 2+ ) arranged ran- 
domly on a pyrochlore structure - than that anticipated 
from a collective paramagnet. 

Recent x-ray absorption measurements 
on Y2M02O7 [0 show that despite the high degree of 



stoichiometry in this material, there is a small but mea- 
surable amount of disorder in the Mo— Mo bond lengths 
frozen into the system at low temperature. In addition we 
note that the liquid like diffuse neutron scattering data 
observed for CsNiCrFg seem to be best described by lo- 
cally collinear spin configurations Thermal selection 
of a collinear state has been considered as a possibility 
in the Heisenberg pyrochlore antiferromagnet. However, 
it has been shown not occur in the disorder-free fic|| and 
site-diluted cases Q. 

In this paper we study the effect of weak bond dis- 
order on the classical Heisenberg antiferromagnet on the 
pyrochlore lattice. We find that bond disorder provides a 
mechanism for the lifting of the the ground state degen- 
eracy, with the selection of a locally collinear state. This 
kind of effect is known as configurational order by disor- 
der |lq |, since the ground states selected by the pertur- 
bation typically incorporate long-range order. However, 
we demonstrate that configurational order by disorder is 
in turn frustrated on this lattice, since the optimal con- 
figuration of individual tetrahedra are not globally com- 
patible, which leads to the collinear order being only of 
short range. In addition, for finite disorder strengths, the 
ground states even of individual tetrahedra are no longer 
necessarily perfectly collinear. A rugged energy land- 
scape arises from the many ways of finding compromises 
between the competing requirements imposed by the dis- 
order. Our Monte Carlo simulations do indeed find that 
the development of local collinearity is accompanied by 
the onset of glassy behaviour. 

The Heisenberg Hamiltonian for the pure pyrochlore 
antiferromagnet is 



H = J ^2 Si ■ Sj = - \L 



NJ, 



(1) 



where J > is the coupling constant, N the number of 
spins and L p is the total spin on the p th tetrahedron. The 
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sum on < ij > runs over all nearest neighbor bonds. The 
ground state condition L p = leaves one degree of free- 
dom per tetrahedron unconstrained, thus yielding a con- 
tinuously degenerate ground state of extensive dimension 
N/2 and localised zero-energy degrees of freedom which 
fluctuate all the way down to T = [|lo| . 

Thermal order by disorder, in which ordered states 
are selected at small but nonzero temperature due to 
the large entropy of the soft fluctuations around them, 
has already been discussed for this model [|l0|,^6). The 
candidate states for selection are the discrete subset of 
collinear ground states, which have the broken uniax- 
ial symmetry of a nematic and are exponentially numer- 
ous in N. However, localising the system in phase space 
near points of this discrete set of collinear states is pro- 
hibitively costly in entropy, even taking into account the 
soft fluctuations around the collinear states: thermal or- 
der by disorder, unusually for frustrated magnets, is ab- 
sent here. 

We add a small disordered term, Hi, to Hq: 

H i = J >\.,S, ■ S, ; -So < S it j < S , (2) 
<i,j> 

with Si j taken from a continous and rectangular distri- 
bution. We show, in the inset of Fig. 1, the energy per 
spin at low temperatures (T/J =0.0006) for different dis- 
order strengths. As 5q increases, the energy per spin de- 
creases linearly, illustrating that the disorder relieves the 
frustration, allowing lower energy states for specific con- 
figurations. One can see that these should be collinear 
by recognizing that the energy on a given tetrahedron is 
minimized if the spins on the weakest bonds (i.e those 
with Sij most negative) are ferromagnetically aligned, 
while those on the strongest bonds are antiferromagnet- 
ically aligned. The energy of a single tetrahedron can 
therefore be reduced compared to a non-collinear state 
by judicious selection of which spins are up and which 
arc down. In the appendix we give a proof that the low- 
est energy state is indeed a collinear state for a single 
tetrahedron in the limit of weak disorder. 

As a test for the formation of collinear order we have 
calculated the probability density function P(w) for the 
scalar product between neighbouring spin vectors, Wij = 
Si ■ Sj, by Monte Carlo simulation. For collinear states 
with 4 antiferromagnetic and 2 ferromagnetic bonds per 
tetrahedron, P(w) is bimodal with delta functions of 
strength 2/3 and 1/3 at w = —1 and w = 1 respec- 
tively. In Fig. 1 we compare P(w) for a system of 1024 
spins in a cubic cell with periodic boundaries for a sys- 
tem without disorder and for a system with So/ J = 0.1. 
In both cases, the data were collected over 10 6 Monte 
Carlo steps per spin (MCS) at T/J = 0.001 and the ini- 
tial configuration generated by cooling in stages of 10 5 
MCS. In the pure case on reducing the temperature a 
peak develops at w = — 1 indicating the presence of weak 
local antiparallel correlations. In the presence of bond 
disorder the low temperature feature at w = — 1 is much 
more pronounced. In addition a clear peak in probility 
is present also at w = 1, as expected from the configura- 
tional order by disorder. Nonetheless, collinear ordering 
clearly remains incomplete. 
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FIG. 1. Probability density function P(w) for 1024 spins. 
The full line is for 6o = 0.1 and the dotted line for So = 0. 
The inset shows the energy per spin with bond energy taken 
from the distribution (Eq. 2). Values are from N = 1024 spin 
systems averaged over 10 disorder configurations. 

To see why this is the case, note that global collinear 
order by disorder is frustrated by the presence of closed 
loops of bonds on the lattice. If the single-tetrahedron 
ground states imply an odd number of antiparallel bonds 
around any closed loop, a global ground state obviously 
cannot be constructed consistently out of those for single 
tetrahedra. 

Collinear order could be retained at the cost of having 
higher-energy configurations on some of the tetrahedra, 
but this strategy is not promising due to the large number 
of closed loops of bonds. The alternative is to introduce 
canting, and this is indeed what happens, as borne out 
by our simulations results. The large number of possi- 
ble compromises imposed by the loops leads to a large 
number of states with similar energy but separated by 
barriers - a rugged energy landscape. 

To leading order in So the local energy minima lie on 
the original ground state manifold with L p = 0, while 
to higher order they move off this surface towards the 
minima of H±. For collinear states the spectrum of exci- 
tations of Hq perpendicular to this surface contains soft 
modes, and therefore canting away from collinearity due 
to higher order effects can be considerable, even for mod- 
est values of So ■ The higher order effects therefore provide 
a mechanism for the destruction of collinearity, even for 
a single tetrahedron and even for modest values of So. 

We have tested for both spin glass behaviour and global 
collinearity, since the ground state manifold is now dis- 
connected due to the formation of energy barriers which 
could induce glassy behaviour and loss of ergodicity. 

The nematic order parameter Q, is the largest eigen- 
value of the traceless matrix [lS| ] 

i 

where a = x,y,z and S a p is the Kronecker delta. The 
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Edwards- Anderson spin glass order parameter |l9| is 

4E<^> 2 - ( 4 ) 
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FIG. 2. Nematic order parameter Q (filled symbols)and 
spin glass order parameters q EA (open symbols) against tem- 
perature, calculated for 128 spins (circles), 432 spins (trian- 
gles) and 1024 spins (squares). 

In Fig. 2 we show Q and qEA against temperature 
for S = 0.1 for N = 128,432 and 1024 spins. The 
data is averaged over 10 — 30 configurations of disor- 
der, with at each temperature 10 5 MCS for equilibration 
and 10 5 for production. The data clearly show the on- 
set of glassy behaviour at a temperature of the order of 
T/J = 5q/J = Q.1. There is a marked decrease of qea{T) 
for T < 0.08 from N — 432 to TV = 1024, suggesting that 
the freezing occurs over microscopic clusters only. More- 
over, preliminary results for the susceptibility suggest, 
at least for small systems, that there is no bifurcation 
between field-cooled and zero field cooled measurements. 
However, the time scales involved do increase by orders 
of magnitude in the presence of bond disorder: the relax- 
ation time at the lowest temperature studied, T/J= 0.02 
is in excess of 10 6 MCS, while for the pure system it is 
of order 10 3 MCS. 

The nematic order parameter scales to zero at low tem- 
perature roughly as Q(T = 0) ~ A/yN confirming that 
the collinearity is only a local phenomena, with correla- 
tions over about 4 tetrahedra. 

We conclude that, although the ground state for a sin- 
gle tetrahedron is perfectly ordered, the state attained 
for the larger system is only partially ordered and is dy- 
namically frozen on the time scales we have considered. 
This may not only be because the true ground state is not 
perfectly collinear, as discussed above, but also because 
a loss of ergodicity in the simulations can mean that the 
system gets stuck in a metastable state. 

Our results therefore provide a plausible mechanism 
for the local collinearity proposed to fit the neutron scat- 
tering data of Harris et al on CsNiCrF 6 || . Bond dis- 
order alone with isotropic Heisenberg spins is probably 



not sufficient to explain the glassy behaviour observed in 
Y2M02O7. However it is sufficient to remove the anoma- 
lous liquid like behaviour of the classical Heisenberg an- 
tiferromagnet and ensure dynamical freezing at a well 
defined temperature, for any finite system. One might 
naively expect that, due to absence of long-range nematic 
order, the present system renormalizes to an effective 
randomly frustrated isotropic Heisenberg spin glass, at 
length scales beyond the short-range correlated collinear 
regions. It is generally believed that the isotropic three 
dimensional Heisenberg spin glass has a zero temperature 
spin-glass freezing transition pp| , with no field-cool/zero 
field-cool split in the magnetization as is also indicated 
by our simulations. Within this picture, the observed 
spin glass transition in Y2M02O7 Q must come from 
anisotropic perturbations Efl. Site disorder has also 
been considered previously |14||, but it does not lead to 
glassy behaviour: substitution of a magnetic ion does not 
destroy the continuous ground state manifold but only re- 
duces its dimension. Bond disorder, on the other hand, 
should lead to energy barriers of a height set by the dis- 
order strength, which is just what we find in Fig. 2. More 
work is clearly required here. In particular, further sim- 
ulations are needed to investigate the nature of the glass 
transition in more detail. It would also be of interest to 
re-examine neutron data with a view to fitting S(q, lo) 
using local collinear spin configurations. 

Finally we note that the longest standing theoretical 
challenge in this field is motivated by the most studied ex- 
perimental system SrCr^Ga^-zOig (SCGO) which con- 
sists of bi-layers of two kagome lattices of magnetic ions 
which are coupled, forming a slice of pyrochlore lattice. 
There is experimental evidence that, at low temperature 
the spins form into a glassy co-planar phase |^,^2|,^3| but 
no theoretical explanation for this result exists. The ex- 
perimental compound is non-stoichiometric and so sub- 
stitutional disorder is inevitable. The thermal order by 
disorder seen in the classical kagome [|ll| antiferromag- 
net has been shown to be extremely sensitive to dis- 
order |Q, with 2% dilution being sufficient to destroy 
coplanarity, and so cannot be responsible for the ob- 
served coplanarity. Further to this, in presence of the 
12fc — 2a — 12k kagome— triangular— kagome coupling, 
the coplanar thermal order-by-disorder found for a single 
kagome layer is destroyed, with or without site or bond 
disorder. It is therefore tempting to speculate that a vari- 
ant of the mechanism presented here could lead to con- 
figurational coplanar order by disorder in parallel with 
glassy behaviour, in agreement with experiment at long 
last. However, we have not, as yet, managed to isolate 
such a variant. 



3 



ACKNOWLEDGMENTS 



This reserach has been funded by NSERC of Canada 
and the Pole Scientifique de Modelisation Numerique 
(PSMN) of the Ecole Normale Superieure de Lyon. M.G. 
and P.H. are grateful to the Association of Universities 
and Colleges of Canada for a travel grant from the Going 
Global STEP program. M.G. acknowledges the Research 
Corporation for a Research Innovation Award and a Cot- 
trell Scholar Award, and the Province of Ontario for a 
Premier Research Excellence Award. 



APPENDIX A: 

We consider the Hamiltonian H = H n + Hi for a sin- 
gle tetrahedron, and treat Hi in first order of pertur- 
bation theory. The ground states of Hq have total spin 
L = 0. They can be parametrized as Si = (0, 0, 1); S 2 = 
(0, sin a, cos a); S3 = (sin/3 sm<f>, sin/3 cos <f>, cos/3); 5*4 = 
— «Si — S2 — S3. The normalisation of S4 imposes the re- 
striction cos a + cos /3 < 0, as well as a condition on </> [[l7| . 
Therefore, 

S\ • S2 — cos a , 
Si - S3 = cos /3 , 
Si ■ S4 = —1 — cos a — cos (3 . 

In a ground state of Hq, the sum of any pair of spins is 
opposite to the sum of the other pair of spins, and hence 

S3 ■ 04 = Si ■ S% , S2 ■ S4 — Si ■ S3 , S2 ■ S3 — Si ■ S4 

The dot products between two spins can thus be ex- 
pressed using cos a and cos/3 only, this means that, re- 
stricted to the ground states of Hq, Hi must have the 
following form: 

Hi/J = Acosa + Bcosf3 + C , (Al) 

where A, B and C depend on the 5: A = 5i 2+^3, 4 — ^1,4 — 
£2,3, B = <S 1)3 + 8 2 a - t>i,4 - 8 2 ,3, and C = -£ 1)4 - £ 2 ,3- 

We now have to distinguish the following cases: (I) 
A 7^ B. Without loss of generality, let A > B. (a) For 
A > 0, B 0, the global minimum of Hi is given by 
cosa = —1 and cos/3 = — sign(S). (b) If A > 0,B = 
0, then cos/3 is undetermined, (c) For A < 0, Hi is 
minimized by cosa = — l,cos/3 = 1. (II) A = B (a) 
A > 0: Hi is minimized by cosa = cos/3 = —1. (b) 
A < 0: Hi is minimized by cosa + cos/3 = 0. (c) A = 0: 
Hi is independent of a, /3. 

From this, we see that two things can happen. Ei- 
ther all spins are collinear, in the cases Ia,Ic,IIa, or the 
degeneracy is not fully lifted, in the special (finetuned) 
cases lb, lib, He. For the latter, a collinear state still ex- 
ists which minimizes the energy. In other words, a non- 
collinear state is never selected by the disorder as unique 
ground state, and the generic case A ^ B, AB ^ has a 
unique ground state which is collinear. 
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